LARGE DEVIATIONS OF THE CURRENT IN STOCHASTIC 
COLLISIONAL DYNAMICS. 
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r-^ Abstract. We consider a class of deterministic local collisional dynamics, showing how 

,_^ to approximate them by means of stochastic models and then studying the fluctuations 

f*-) of the current of energy. We show first that the variance of the time-integrated current 

£SJ is finite and related to the conductivity by the Green-Kubo relation. Next we show that 

,__! the law of the empirical average current satisfies a large deviations principle and compute 

explicitly the rate functional in a suitable scaling limit. We observe that this functional 

is not strictly convex. 
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1. Introduction 



The paper J18] introduced a class of stochastic dynamics aimed at modeling Hamiltonian 
dynamics describing aerogels. Those stochastic dynamics describe also pretty accurately 
the behaviour of deterministic models made of tracer particles and fixed scatterers [11| IT^j 

CSJ El SHI EI] • In non-equilibrium statistical mechanics, it is now a well established fact that 

the large deviations functional of the currents of the relevant conserved quantities play a 
role similar to the thermodynamic potentials in equilibrium, see [U E] for general overviews 

\^q of the subject. In this paper, our goal is to go one step further than in [TB] and study 

the fluctuations and large deviations properties of the energy current of one of the models 

[^ considered there: the confined tracers. A confined tracer is a particle that moves freely in an 

interval and is reflected at the boundaries of the interval with a random speed v distributed 
according to the distribution: 



„2 



•** ipp(v) = five P v * , v > 0. 

rS 

This random reflection models the action of a "hot" wall and /3 is the inverse of the temper- 
ature of the wall where the collision has taken place. We first explain the formal connection 
between the deterministic and the stochastic dynamics. In particular we derive the ex- 
pression of the infinitesimal generator of the stochastic dynamics, that is an example of 
piecewise deterministic dynamics. Next, the Green-Kubo relation for the stochastic dy- 
namics is rigorously derived. 

Finally, we obtain our main result: we show that the law of the empirical average current 
satisfies a large deviations principle and compute its scaling limit. A striking feature of the 
rate function is that it contains linear (affine) parts. Namely, we obtain the following limit 



Date: July 22, 2010. 

Key words and phrases, collisional dynamics,heat conduction, large deviations of the current, scaling 
limit, renewal processes. 

1 



R. LEFEVERE, M. MARIANI, AND L. ZAMBOTTI 



rate function: 



G(j,r,T):={ 



0- KT ) 2 if 7 - T > KT 2 
if JT G [0, KT 2 ] 

£ if jre[ _ Kr 2 i0] 



2T 



(1.1) 



-KT 



. 4kT 2 ir -? T < 

where r is the rescaled temperature gradient applied at the boundaries of the interval and 
T the arithmetic mean of the temperatures. The graph of Q as a function j is given in 
Figure 1 below. In an accompanying paper [19] . we discuss the implications of this result 
for the macroscopic fluctuation theory [21 [3j HI El El [7] of deterministic systems modeling 
aerogels. A similar large deviations functional has been found in the context of random 
walks in random environments [U [15] . 
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Figure 1: Plot of Q as a function of j for kt = kT 2 = 1 



In order to obtain the large deviations principle and the rate function, the main difficulty 
to overcome is the lack of strict convexity of the candidate rate functional G(j, t, T), which 
is obtained as the Legendre transform of a cumulant generating function. In this situation 
the Gartner-Ellis theorem does not yield the full large deviations principle and we have to 
use a more detailed analysis, in particular in order to obtain the lower bound. 

2. From deterministic to stochastic dynamics. 

Consider N particles of unit mass with positions and momenta (q, p) = {(qi, Pi)}i<-<w> 

with qi,pi G M. d . The positions are measured with respect to N fixed centers located on a 
ID lattice. The Hamiltonian H takes the form, 



N 



ff(p,q) = £ 



PT 

2 



+ V(qi) + U(qi-q t+1 ) 



(2.1) 



where the interaction potential U is equal to zero inside a region Qjj C R with smooth 
boundary A of dimension d—1, and equal to infinity outside. Likewise, the pinning potential 
V is assumed to be zero inside a bounded region £ly and infinity outside, implying that the 
motion of a single particle remains confined for all times. The regions Qjj and fly being 
specified, the dynamics is equivalent to a billiard in high dimension. A typical example 
of the dynamics we wish to consider is given by the figure below. The circles move freely 
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within their square cells and collide with each other when they both get sufficiently close 
to the hole located in the wall separating two adjacent cells. 









Figure 2: Simplified one-dimensional aerogel dynamics 



Physically, those models describe aerogels, i.e. gels whose liquid components have been 
removed and replaced by atoms of gases. In this section, we describe how to replace the 
interaction of one particle with its nearest-neigbors by an interaction with local stochastic 
heat baths. The model is an approximation based on macroscopic fluctuation theory and 
a local equilibrium picture: the action of the whole system on a single atom through its 
nearest neighbors is replaced by the action of infinite thermal baths on which the particle 
motion itself has no direct influence on the microscopic time scale |19j . One should expect 
the approximation to hold true whenever numerics show that the systems are close to local 
equilibrium, i.e. in a weakly interacting regime, see \1'6\ I14|. This approximation should 
also apply to some models considered in [UJ HI E3 1201 EI] as explained in |T9j . 



2.1. Deterministic dynamics. We denote by Ajj and Ay the boundaries of the regions 
ftu et Qy. We define also: 



n = {(qi,...,qjv)| (q» -qi+i) e o 



u, i 



1,...,N-1, qi G Vt 



y, I 



1,...,N}. 



Example. The simplest example of the type of dynamics we are interested in is given by 
a limit of models with smooth interaction potentials. In the definition of the Hamiltonian 
(2.1), take d = 1 and replace the potentials V and U by V& and U^ where V^{x) = /&(§) 
and Uk(x) = /&(£), where fk(x) = x 2k /2k. When k — > oo, the limit interaction potentials 
Uoq and V^o are square wells of infinite heights whose walls are located respectively at ±6 
and ±a. 



V' 



+oo if \x\ > b 
if \x\ < b 



U a 



+oo if |x| > a 
if \x\ < a 



(2.2) 



Each particle on the lattice moves freely on a one-dimensional cell of size 26, changing 
directions at the boundaries. The interaction between a pair of particles acts when the 
difference between the positions of the two particles reaches the value a, at which point 
they exchange their velocities. 



2b 
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Figure 3: The complete exchange model 
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The motion of a given pair of particles at sites i, i + 1 is described as the motion of a point 
particle on a two-dimensional billiard table described by 

Q = {(x\,X2) & R , \x\\ < b, \x2\ < b, \xi — %i\ < a} • 



2b 





2b -a 



Figure 4: Billiard giving the motion of two particles in the complete exchange model 

It is straightforward to build analogous models on higher-dimensional lattices. Those 
models have been introduced in [22] and called the complete exchange models [14J. 



We denote by p = (pi, . . . , pjy) and q = (qi, . . . , q^r) the vectors made of the momenta 
and positions of each particles. The evolution of a probability distribution P t (p,q) over 
phase space is given by 



(2.3) 



dtPt = ~ Eili Pid^Pt, if q G O 
P t (qi,...,qAT,pi,...,pjv) = Oif q ^ Vt 

with specular boundary conditions: 

Pi(q,p)=P t (q,p-2(n-p)n), if q G dQ 

and n is a normal vector to the boundary at point q. We denote by /i(q, p, t), the marginal 
probability distribution of the i-th particle 



fi(q,p,t) = Y[d Pjd qj P t (qi, . . . ,qjv,Pi, ■ • • ,Pn) 



(2.4) 



and similarly /ij(q, p, q',p't) denotes the probability distribution relative to the momenta 
and positions of particles i and j. The equation for the evolution of the set of probability 
density of a single particle is given in terms of collision with the boundaries of the walls of 
its own cell and with its neighbors 



d 

~di 



/i(q, P, t) = -p • V q /i + L w /i + L c (fu +1 ) + L c (f,, 



i.i— 1 , 



(2.5) 



Here L w accounts for the collisions of the particles with the walls of their respective cells: 
L w (/ 4 )(q,P,) = <Mq)( P • n) + [/ i (q,p - 2(p • n)n) - /,(q,p)], (2.6) 

here and below n denote a generic unit normal vector to the boundary Ajj or Ay at point 
q. We use the notation 

„± _ / l x l ^ ± x > 
if ±x<0. 
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The collision term L c (fa±i) for the collisions of the i-th particle with the i ± 1th is given 
by, 

L c {fi,i±i) = / dp a dq / 5 A (q-q / )(P" L -pf) + [/i,i±i(q,Pc,q / ,Pb)-/i,i±i(q,P,q / ,Pa)] (2.7) 

with p^ = p x , Pc = Pa '> Pc — Pc*i = P — p n, and pt, — p^n = p a — p^ n. One can check 
that the distribution 

N 

P cq ( qi , . . . , dN, pi, . . . , p N ) = Z-^nCoi, , Qiv) J] e_/3p?/2 ( 2 - 8 ) 

8=1 

is stationary for any inverse temperature j3. We will explain at the end of the next subsection 
how to thermalize those systems at different temperatures at their boundaries. 

Example. Coming back to the example of the complete exchange models, we see that the 
different collision terms take the form 

L w (fi)(q, P ) = 5 ±b (q)p ± [f i (q, -p) - f t (q,p)}, (2.9) 

L c {fi,i±i)= fd P 'dq'd ±a (q-q / )( P -p') ± [f iA±1 (q,p\q > , P )-f i , i±1 (q,p,q > , P % (2.10) 

As our goal is to see how to approximate those dynamics by stochastic ones. We now 
describe in detail the case of a free particle confined between two hot walls. The action of 
the walls is purely random and express the effect of the contact of the particle with a very 
large system whose temperature is fixed and constant. 

2.2. A free particle confined between two thermal walls. Let us consider a particle 
moving in the interval [0, 1] and reflected at the left and right boundaries boundaries with a 
random velocity p whose absolute value is distributed according to distributions respectively 
(j)p L and (f>p R , with: 

<f>p(dp) = f3pe-^l {p>0) dp. (2.11) 

The temperature of the right, respectively left, wall is Tr := 1//3r, resp. Tj_, := l//3^. As 
the particle is reflected at the boundary of [0,1], we understand that the distribution of 
velocities at the thermal walls has a sign which is opposite to the sign of the velocity of the 
incoming particle. 

In the following, we recall some notation and results following |18| . Let E = {— 1,+1}, 
(qo,Po) the initial data of the particles, and define o"o = sign(po), <?k = {— l) fc co for k > 0. 
For a £ E. In |18j a more general case has been addresed, where (a^) is Markov chain on 
the state space E. 

As explained in [18], this is a Markov process, the notations are similar to [18!. We 
explain here the special case corresponding to the case we are interested in. The state 
space of the associated Markov chain is E = {— 1,+1}. Let (qo,Po) the initial data and 
velocity of the particle. We define o"o = sign(po)- We consider now the Markov chain 
(cfe)fc>o i n E with initial state Xq = <jq. In fact, the Markov chain has a deterministic 
evolution cr^ = (— l) fc o"o, k > 0. 

For each a G E, we write a = ^ (c + 1)- Then the time of the first collision with a wall is 

So = S (q ,p ) := > 0, 

Po 
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We now define the sequence of times the particle takes between two subsequent visits to the 
scatterers. Conditionally on (o"fe)fe>o, the sequence (rk)k>i is independent with distribution 
defined by 

Pcfc-l f Po h _ x 



(r k G dr | a k -i) = -^- exp I — °^ ) 1 (r>o) dr, (2.12) 



where /3_i = /Sz, and /3i = (3r. In other words, the conditional law of l/r k is (j)g a : 

P(l/r fe G dp | <Tife_i) = / 3 eTfe _ 1 p e-^-i ^ l (p>0 ) dp. (2.13) 

The time of the (A; + l)-st collision with one of the two walls is 

Sk:=S + Tl + --- + T k , k>l. (2.14) 

Notice that (S k ) k >o is a standard renewal process, see [T|. Before time So, the particle 
moves with uniform velocity po- Between time Sk— i and time S^, the particle moves with 
uniform velocity ^ and (S k )k>o is the sequence of times when q t G {0, 1}. In particular we 
define the sequence of incoming velocities Vk at time Sk 

v :=p , v k := — , k>l. (2-15) 

Tk 

We define the stochastic process (qt,Pt)t>o with values in [0, 1] x R* 

{qo+Pot,Po) if t < S , 

(2.16) 



(qt,pt) ■-- 



a k -i + f (i - 5 fc _i), a ) if S k -x <t<S k , k> 1, 



The energy exchanged between the two walls during a time interval [0, t] is given by 



J[0,t}:=l Yl 



2 ^ V k Gk - 
k>l:S k <t 



and we have shown |18] that 



iimij[o,t] = K(r L -r fl ) (2.17) 

4— >oo I 

where Tr = 1/Pr, Tl = 1/(3l and k is by definition the conductivity of the model, given 
by: 

We will establish a correspondence between the deterministic dynamics and stochastic dy- 
namics by using simplification of the evolution equation of the probability distributions. 
Thus, we need to write the infinitesimal evolution of probability distributions under the 
stochastic dynamics described here. We set for bounded Borel / : [0, l[xM+ i— >• R 

Ptf(qo,Po) :=E(f(q t ,p t ))=E(f(F(q,p,t,(T n ) n > 1 )), (q ,p ) G [0,1] x R* . 

In the appendix, we prove the following: 



LARGE DEVIATIONS OF THE CURRENT IN STOCHASTIC COLLISIONAL DYNAMICS. 



Proposition 2.1. For all f,g : [0,1] x R + 

derivatives: 



bounded with bounded continuous first 



d 
dt 



I dq I dpg(q,p)P t f(q,p) = 

dqg(q,p)pf q (q,p)+ dppg(l,p) 4>p R {du) (/(l, -u) - f(l,p)) 



dp 



d P pg(0,p) / <t>p L (du)(f(0,u)-f(0,p)). 



+ 



From this computation we obtain a formal expression for the generator: 



Of 

Lf{q,p) =P~q-+P S {dq) 

+ p + 5 1 {dq) 



And similarly, one may obtain an expression for the formal adjoint: 

dg 



[ ^ L (du)(f(0,u)-f(0,p)) 

JR+ 

cf>p R (du)(f(l,-u)-f(l,p)) 



L*g(q,p) = -p 



0q 



o~o(dq) 
Si(dq) 



P g{Q,p)-<t>p L (p) I duug(0,u) 
P~ g{hp) - 4>]s R (p) / duug(l,u) 



(2.19) 



This tells us also how to describe the thermalization on the boundaries of the deterministic 
dynamics described in the previous subsection. For instance, in the case of a system of N 



particles described by the complete exchange dynamics described in (2.9) and (2.10), one 
simply replaces L w (/i) and L w (/n) by: 



and 



L h {h){q,p) = -5± b {dq) 



L l5R (f N )(q,p) = -S± b (dq) 



P fi(q,P)-<Pp L (p) duufi(q,u) 



P ± fN(q,p) ~ 4>%{P) I duuf N (q,u) 



(2.20) 



(2.21) 



where (pg(p) = /3(p) =t e - ^^ . In more general models, the thermalization is similar: the 
specular reflections of the particle on the walls of its own cell is replaced by the action of a 
thermal wall. 

We notice here that our process is piecewise deterministic, at least in the sense that 
randomness acts in a discrete (random) set of times and the motion is deterministic in 
between. Piecewise deterministic processes have been extensively studied, see for instance 
[16], however our case does not fit in the standard framework. Indeed, in the literature 
one finds piecewise deterministic processes with generators like (2.20) and (2.21) with the 



measures S±b(dq) replaced by some function on the state space; this corresponds to a noise 
which can act at a positive and finite rate all over the state space, while in our situation the 
noise act every time that, and only when, the process hits the two lines {{q,p) : q = ±1}- 
Therefore the standard theory can not be applied to our processes (2.16). 



R. LEFEVERE, M. MARIANI, AND L. ZAMBOTTI 



2.3. Stochastic approximation of the deterministic dynamics. We want to describe 
the dynamics of a given particle when its neighbors have positions and velocities distributed 



according to a local equilibrium distribution. In order to do so, in the expression (2.10), 
we set: 



fi,i±i(q,p,q',p') = fi(q,p) 



6 



i±\ 



8b 2 n 



-ft±i 



"k[-b,b](d)- 



(2.22) 



Then (2.10) becomes 



Pi±i 
8b 2 TT 



L c (f i , i± i)(q,p)=\l'^ : t hb , b - a] (q) dp\p-p'y[f l (p',q)e-^^ -f l (p,q)e-^ P -] 



+ l[«-6,6](g) dp'(p-p')-[fi(p,q)e-^ ±1 ^-fi(p,q)e-^ 



(2.23) 

We recall that b < a < lb. Heuristically, the interpretation of the dynamics with the 
collision term modified as above is easy to describe: each particle moves freely and bounces 
back and forth in its cell of size 26 except when it enters two "critical" regions of size 
2b — a located near the boundaries of the interval [—6, b]. There, at a random point it may 
collide elastically with a particle whose velocity is Maxwellian. Simplifying further, one 
may contract the critical regions to the points at the boundaries of the interval and replace 
the collision with a Maxwellian particle by the collision with a thermal wall. One is then 



lead to equations of the form (2.20), (2.21). One may of course use the same strategy for 
any of the deterministic collisional dynamics we described above. The geometry may be 
completely different but the idea is always similar: each particle moves freely in its cell 
except in a "critical region" where collisions with neighbors may occur. If one is interested 
in the local equilibrium dynamics, one is led to a dynamics described by free motion and 
interactions with thermal walls. Pictorially, the model of figure 2 is transformed into the 
model of figure 5 below. 




Figure 5: Collisional regions have been replaced by interactions with thermal walls 

Clearly, the direction orthogonal to the array of heat baths is irrelevant to the transport 
of energy and one is thus lead to the model of confined tracers introduced in [18J. In other 
words, one replaces the effect of the neighbors of a given particle by the effect of heat baths. 
The dynamics of each particle is given by the motion of a particle confined between two 
thermal walls as in subsection 2.2 We discuss now the large deviations properties of the 



current carried by a single particle confined between two hot walls. 



3. Fluctuations and large deviations of the current. 

Our final goal is to study a suitable scaling limit of the large deviations functional of the 
energy current in the case of a confined tracer. In order to identify the different physical 
quantities appearing in the scaling limit, we start off by computing the variance of the 
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current in equilibrium 



lim V q ((J[0,t]) 2 ) 



where the expectation is taken with respect to the process defined in subsection |2.2| with 
Pl = Pr = P = T^ 1 for some T > 0. In principle this quantity should be computable as 
the second derivative with respect to A of the cumulant generating function 

/(A):= lim JlogE cq (exp(AJ[0,t])), (3.1) 

Unfortunately, this amounts to assuming that one can exchange the limit t — > oo and the 
derivatives with respect to A. As we are not aware of any general argument justifying this 
exchange we compute explicitly the variance of the current below. We find that it coincides 
indeed with the second derivative of the cumulant generating function obtained in [18J. We 



will see that the variance (3.3) appears in (4.2) 



3.1. Variance of the current in equilibrium. We recall that we can associate to the 
renewal process (Sk)k>o the counting process 

oo 

fc=l 
Proposition 3.1. The Green-Kubo relation 

<r(T) := lim V q ((J[0,t]) 2 ) = 4k(T)T 2 (3.3) 

t— >oo t 



holds true, where k(T) is given by (2.18) when P L = P R 



Proof. We recall that the energy exchanged between the two walls during a time interval 
[0, t] is given by 



1 1 Nt 

J[o,t] = - E *o(-i)*«2 = «I>(-i) 



Write 



Then 



2 ^-^ 2 

k>l:S k <t k=l 



HJ[o,t]) 2 = (f^(-i) k vA 



N t 

v k 



/N t N t 

4E cq ((J[0,t]) 2 ) = E cq £(-l)*t£ X>l) fe vl 

\k=l fc'=l 



^fc+fc' 2 2 

v k v k> 



E E Eeq (l(JV t =0 (-1) 

£=1 l<k,k'<e 

oo e 

EEMV*=W)+ 2 E E Eeq(v«=o(- 1 ) 

£=1 fc=l £=2 l<fc<fc'<^ 



v k v k' 



Observe now that conditionally on (Nt = £), the variables {vi,k < £} are exchangeable. 
Therefore 

oo I oo 

E E E ^ ( V«=o O = E ^ E ^ I 1 w=o u i) = E ^ (^ u i) • 

£=1 fe=l £=1 
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Moreover, by exchangeability 



£ £ E eq (l (JVt= ,) (-!)*+*'« 

^=2 l<fc<fc'<€ 



X' 



E E (-l) W 'E eq (l(iv t =^N2 2 ) 

1=2 l<k<k'<£ 



^(-l)L^J E cq (l (iVt= ^ 2 2 ) = E cq ((-l)WJ vfv 



and therefore 

4E cq (( J[0, t]) 2 ) = E cq (N t vf + 2 (-1)L"*/2J y2 ^ 

We now compute 

lim ^E cq ((J[0,t]) 2 ) = lim ^(N t u\ + 2{-l)^ N ^vlvl 

where 

// = /3 / e - ^ (if 



(3.4) 



//3vr 

JO v T 

and the last equality follows because the second term is uniformly bounded in t and because 

lim E cq ( — - - = 0. (3.5) 

In order to see this last point, it is enough to prove that Var(-p) — > 0, since we already now 
that Nt/t — > 1//J, in L l by the Renewal theorem. We approximate the variables {jk)k^^ by 
a sequence of truncated variables (t^^n, r^ = r^ A n. Each r^ - has a finite variance, say 
a\ and average [i n and theorem V.6.3 in [1] implies for the corresponding renewal process 
N? that 

lim -Var (NJ 1 ) = % 

t^+oo t Hi 

and thus for each n, 

/ /v ra i \ 2 
lim E eq -«- = 0. 

Notice now that tJ} < r^ a.s. for all n, fc implies iVf < iV t n a.s. and therefore 

Var (JV t ) = E ((N t ) 2 ) - (E (iV t )) 2 < Var (JVf) + (E (iV^)) 2 - (E (iV 4 )) 2 . 
Then for all n > 1 

lim -jVar(iV t ) < -=- - -^. 

Since // n = E(r{ 1 ) f E(n) = \i by monotone convergence, then we obtain that indeed 

Nt 

t 



Var(-p) — y and we have proven (3.5). Now 



r'OO 2 Q 

E(v$) = P I v 5 e~^dv=^. 



LARGE DEVIATIONS OF THE CURRENT IN STOCHASTIC COLLISIONAL DYNAMICS. 11 

Therefore 

lim V q ((J[0,i]) 2 ) = * x [f = AT 2 K (T). (3.6) 

□ 



3.2. Cumulant generating function. We state first the results of [18] which apply to 
the case of a single tracer confined between to thermal walls as described in subsection |2.2| 
We define B := (/3_l,/3r) and 

f(X,P L ,P R ) = f(X,B):= lim J logE(exp (AJ[0,t])) , A G ] - (3 R , (3 L [, (3.7) 

r+oo 2 

C(x,rj):= ve~«~ x ^dv, v,x>0, (3.8) 

Jo 

and 

F(\,r ] ,B):=(3 L p R C(P R + \,r ] )C((3 L -\, V ). (3.9) 

From [IS] , we have the 

Proposition 3.2. Suppose that (3l < /?_R- 27*e function f(-,B) is convex and continuous 
over } — /3r, /3l[ and satisfies the Gallavotti- Cohen symmetry relation 

f(\,B) = f(p L -p R -\,B). (3.10) 

f(-,B) is analytic on ] — /3 r ,/3l — /5r[U]0,/3l[. Moreover 

(1) V A G ] — (3 R , [3l — f3 R [U]Q, (3l[, f(X,B) is given by the unique solution r/o > to the 
equation 

F(X, m ,B) = l. (3.11) 

(2) //AG [/3 L - fo, 0], ften /(A, B) = 0. 

(3) If\$]-p R ,p L [thenf(\ t B) = +oo. 

(4) §{(0+,B) = K (T L -T«) 

(5) ^ :]0,(3 l [^]k(T l -T r ), +oo[ and ^ :]- (3 R ,f3 L - f3 R [^]- oo,-4T L -T R )[ are 
bijections. 

Proof. The result is proved in |18|. Propositions 4.5, 4.7] (notice however that there is a 
difference of sign, i.e. in |18j one finds statements and proofs are given for /(— A, B)). It 



only remains to prove point (5). We consider A g]0,/3l[. Then by (3.10), f n (X,B) > is 



defined by the relation F(X, f n (X,B),B) = 1. Hence by the implicit function Theorem 

^(X,B) = ~(X,f(X,B),B)(^(X,f(X,B),B) X 

A computation yields 

OF 

^(A,,,B) = 

= PrPl [ v lV2 l {vl - v\) exp (-H - H - tf R + \) V ±-(I3 L - X) ^) dv x dv 2 , 
Jm.1 2 V v i v 2 2 2 ) 
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and 

dF 

= PrPl [ (vi + v 2 ) exp ( -1 - 1 - {p R + \) v l-(i3 L -X) v ^) dvx dv 2 . 

JWLl \ V! V 2 2 2 

Since f(0,B) = 0, by letting A | we find 



dX 



(0 + ,B) 



2 2 

PrPl J R 2 ui ^2 2 («2 - v i) e-Px^t-^^ dvx dv-i 



l l 

fa Pr 



T.ie{L,H] A / K+ exp {-PiT) dv 
--k{T l -T r ). 



— + \l — 

On the other hand, as A f Pl we find with a change of variable that 

— (\,f(\,B),B)~C 1 ((3 L -\y 2 , —(\J(\,B),B)~C 2 ((3 L -\)- 1 , 
for some positive constants C\,C 2 , and therefore 

In order to prove that ^ : ] — Pr, Pl — Pr{^~ >} — °o, —k(Tl — Tr)[ is a bijection, it is enough 



to apply the Gallavotti-Cohen symmetry relation (3.10). □ 



The Legendre transform of the cumulant generating function. We define the Legendre trans- 
form of f(-,B) 

I(j,B)=sup{Xj-f(X,B)}, jeR. (3.12) 

A 

We have the following simple 

Lemma 3.3. The function I(-,B) is positive, convex and finite on K. Setting j* := k(Tl — 
Tr), 

(1) I(-,B) is smooth and strictly convex overM\ [—j*,j*] 

(2) I(j,B)=0forallje[Q,f} 

(3) I(j,B) = (Pr - p L )\j\ for all j G [-j*,0]. 

Finally, I(-,B) satisfies the Gallavotti-Cohen symmetry relation 

I(j, B) = (p L - p R )j + I(-j, B), je R. (3.13) 

Qualitatively, I(-,B) has the profile pictured in Figure 1. 

Proof. By Proposition Q5), for all j € R \ [-f,f] there exists A £] - Pr,Pl[ with 
§{(A,£) = j. By HQ1 Lemma 2.3.9(b)], for such (j, A) we have I(j,B) = jX - f(X,B) and 
moreover A is an exposing hyperplane for j, yielding smoothness and strict convexity. 

Let us now fix j G [—j*,j*] and consider the function g(X) = jX — f(X,B), X G R. By 
Proposition's), g'(X) > for A G ] - Pr, p L - p R [ and g'(X) < for A G ]0, p L [, so that 
the supremum of g over M. is the same as the supremum of g over [Pi — Pr,0]. But on 
[Pl ~ Pr,0] we have f(-,B) = 0, so that g(X) = jX. Therefore I(j,B) = g(0) = if j > 
and I(j,B) = g(P L - p R ) = (Pr - p L )\j\ if j < 0. 
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The Gallavotti-Cohen symmetry relation (3.10) for / becomes now ( |3.13 ), since 
I(j, B) = sup{Aj - /(A, B)} = sup{Aj - f((3 L -f3 R -X, B)} 

X X 

= sup{(/3 L - i3 R - X')j - /(A', B)} = (& L - R )j + I(-j, B). 

X' 



□ 



3.3. Large Deviations of the current. The main result of this section is the following 

Theorem 3.4. Let T > and r > 0. The law fit of J[0,t]/t satisfies a large deviations 
principle as t — )• +oo with speed t and rate I(j) := I(j,B), j G R, i.e. for any Borel set 
A CM. we have the upper bound 



and the lower bound 



lim - log fit(A) < — inf I, 

t^+co t A 



lim - log ^(^4) > —inf/. 

t— >+oo t A 



(3.14) 



(3.15) 



Proof. With no loss of generality let l3l < Pr- Note that I(-,B) is a continuous, convex, 
coercive functions. Therefore the upper bound (3.14) is a direct consequence of (3.7), (3.13) 



and the Gartner- Ellis theorem, see |X0(, Theorem 2.3.6(a)]. 

We first prove that the lower bound (3.15) holds for all A C R \ [— j*,j*], where, as in 

pis] ), 

Pi' ~ PR' 



k(T l - T R ) 



ttJ3l 



+ 



rc/3 R \ - 



Indeed, still by the Gartner-Ellis theorem, see |101 Theorem 2.3.6(b)], the lower bound 
holds for open subsets O C R suc h that, for all j £ O there exists Xj g] — (3r,Pl[ with 
g^ (Xj,B) = j. By Proposition 3.2-(5), this holds true for all j G R\ [—j*,j*] and thus the 
lower bound holds for open (or, equivalently, Borel) sets Ocl\ [~j*,j*]- 

In order to prove the full lower bound for all open subsets of M, we will show the following. 
For each j G [—j*,j*], there exists a sequence (ytjt of probability measures on R such that 



vt — *■ 6j and 



lim sup - H(i/ t | fit) < I (J). 

t— >+oo t 



(3.16) 



Indeed, if (3.16) is proved, then we argue as follows. Let j G [— j*,j*] and let O be an open 
neighborhood of j. Then, using log ft (O) < and Jensen inequality 

logMO) = log f pdv t = log(^- f p-dvt) +logv t (0) 
Jo dvt \M°) Jo dut J 

dfi t 



> log 

Since xlogx > — e _1 for x > 

log MO) > 
> 



i 



v t{0) Jo dv, 



dvi 



-^)/o 10g (S.''"'' 



1 



MO) 
i 

MP) 



H (^t \ih)+ log ( — ^ ) — ^ dfi, 
J c \dfitj dfi t 



H(z/ t \nt)-e 



-i\ 
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Since j G O, v t — *• <5j and O is open, then z^(O) — >• 1 as £ — )• +00. We obtain by (3.16) 



lim - log nt(0) > - lim - R(v t \ /j, t ) > -I(j). 
Therefore, optimizing over j G [—j*,j*] 

lim - log nt(0) > - inf I{j). 

t^+00 t jeOn[-j*,j*] 

Finally, since we know that the lower bound holds on open subsets of R \ [—j*,j*], for a 
generic OcR 

lim - log Ut(0) > max ( — inf I(j),— inf I(j) ) = — inf I(j). 
t->+oo t V 3eO\[-j*,j*\ jeOn[-j*j*} J jeO 



Proof of (3.16) for j G [0,j*].Let us denote 4> + := </>p L , (j) ■= (j>p R . Let us set a := 
j/j* G [0, 1] and for a G { + , -} 

and for e > 

p(p) - p (p) := (1 - a) ^_ + a ^J-^y (3-17) 

and 

^ CT (cip) := p(p) 1 a (dp), r := -J—p^(d P ). (3.18) 

Notice that 7r CT — v a^ a + (1 — a)5o while 7f cr — *• (j)° , as e 4- 0. In particular, the weak limit 
of -K a as e i depends on a, while the limit of 7r CT does not. This property is crucial to the 
argument, in particular in (|3.21) below. 



We denote by P the law on (M+ ) such that under P the sequence (tj)j>i is independent 
and 

(1) for all i G 2N, 1/r, has law (\T 

(2) for all i G 2N + 1, l/n has law </>+. 

and by P^ the law on (M 1 ^*) such that under P^- the sequence (tj)j>i is independent and 

(1) for all i G 2N, 1/t% has law 7r~ 

(2) for alH G 2N + 1, l/r* has law n+. 

For 77 > we also define T t := L - + K 1 ^ 7 ! (p) t j ; and let P 4 '^ be the law on (Rf ) such 
that under P*- 6 ^ the sequence (tj)j>i is independent and 

(1) for all i = 1, ... ,Tt and i G 2N, 1/Vj has law tt - 

(2) for all i = 1, . . . , T t and i G 2N + 1, 1/r, has law vr + 

(3) if i > T + 1 and i G 2N, then 1/tj has distribution 0~ 

(4) if i > T t + 1 and i G 2N + 1, then 1/r, has distribution </>+. 

Let us denote by vt,e,r\ the law of J[0, t]/t under P*' 6 ^. Let us prove now that 

lim lim v* e „ = 5,- . (3.19) 

40 tt+oo ' " J 

By the law of large numbers, under P^- we have a.s. 

,. S Tt .. S Tt T t 7r+(p) + vr-(p) 2(1 + 77) ! , ^ 1 

hm = lim — = — -— = 1 + 77 > 1. 

i->+oo t t^+00 T t t 2 7T+ (p) + TT' (p) 
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However Sx t has the same law under P^ and under W t,e,r) , so we obtain for any e, r\ > 
lim P 4 '^ (S T , < t) = lim P* f — < 1 J = 0. 

t->+oo i->+oo I f J 

or equivalently 

lim ftevm ) = i. 

where 

D t>n := {S Tt > t} . 

Recall {S n > t} = {Nt + 1 < n}, so that on A,r? we have Nt + 1 < Tt ■ Therefore for any 
e, rj > 



lim W*> e > r i(\J[0,t]/t-j\>e)< lim P*({| J[0,t]/t - j\ > e} n A,,) + P t ' e ' r? (A c J 

= lim" Pjr({|J[0,t]/t-j|>e}nA,,) 

t— >+oo 



(3.20) 



By the law of the large numbers and by the renewal theorem [H Proposition V.1.4], we 
have Pjf a.s. as 1 1 +o° 



f 



1 Tf n +i _ ll v\, 2 r n fc+i_ liV i * v\, 2 c nfc+i , 1 e #(tj -^ 

-^0,tj---^^(-l) _- T _^,, 4 ( 1) ->- E,(r 1 + r 2 ) 

_ 1 ^ + (p 2 ) -^~Q? 2 ) _ 1 k + (p 2 ) - n- (p 2 ) 
~ 2n+(l/p)+n-(l/p) ~ 2 ^y + ^y 



and as e 4 

l7f+(p 2 )-7f-(p 2 ) 10 + (P 2 )-^(P 2 ) 



2 J^ + 



7T+(p) 7T-(p) 



2 f, , + ! 



a- 



i fl-i 



0l " /3 



07+ (p) Q7"(p) 



TTfe \ 2 + ^h! 2 



- = af = j, (3.21) 



so that by (3.20) 



lim lim P # (I J[0, t]/t - j\ > e) = 

40 *t+oo 



(3.22) 



which yields (3.19) 



Now we estimate the entropy 
lim" \ R{n, e , v I fit) < lhh \ H (P 4 -^ I P) 

Tt 



lim - J2 (l(ie2N+i) H(fr+ I <f>+) + l (ie2N) H(tt | <fT)) 



(3.23) 



tt+oo £ z 



j=l 



so that 



hm lim - i±(vt, e ,T} Ht) < hm _ - 



(3.24) 



Now, recalling that p := ^qr by (3.17)-(3.18), a tedious explicit computation based on the 
definition of ix a shows that 

H(^+|(/> + ) = / log(-p- J dfr + = _l O g(0(p))+ / logpd# + -> -loga + loga = 
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as e 4 0, and analogously for H(-7f \<p ). Moreover, arguing as above 

2(1 + 77) 



lim lim 

??4.0 e|0 7T+(p) + TT~[p) 



2an < +oo. 



(3.25) 



Thus 



lim lim - H(z/ t £ „ | /z t ) 

ej.0 tt+oo t ' 



which implies, together with (3.22), that there exists a map t h-> (e(t),r/(t)) vanishing as 
£ 1 +°° such that fj := Vt,t(t),r}(t) ~* $j an d nm t ^ _1 H(z^ | /Jt) = 0. 



Proof of (3.16) for j G [— j*,0]. Now the strategy is very similar, but we reverse the 
current constructed for j G [0,j*]. This mirrors the Gallavotti- Cohen symmetry relation 



(3.13). 



g(l+g) 



Set now a := -j/;* G [0, 1], vr ff as in pl7Q-pl8Q and T t := L ^+gffi^ tJ with 7/ > 0. 
Let us define by P i ' <E ' ,? the law on (R5 ) such that under IP*' 6 ' 7 ? the sequence (tj)j>i is 
independent and 

(1) for all i = 1, . . . , T-t and i G 2N, tjj has law tx + 

(2) for all i = 1, . . . , Tf and i G 2N + 1, tjj has law 7r~ 

(3) if i > Tj + 1 and i G 2N, then t> j has distribution <p~ 

(4) if i > Tf +■ 1 and i G 2N + 1, then Vi has distribution + . 



Let us denote by vteri the law of J[0,t]/t under f t ' € ' r i. Arguing as in the proof of (3.19) 
we obtain that under P#, a.s. as i f +°° 



1 



J[0,t] 



1 n~ (p 2 ) - n+ (p 2 ) _ 1 n + (p 2 ) - n~ (p 2 ) 

2 7r+(l/p) + 7r-(l/p) ~ 2 



1 1 1 

7T+(p) 7T-(p) 



and as e i 



l+"(p 2 )-++(p 2 ) 
2 



1 



+ 



1 



-aj 



7T+(p) 7T (p) 



so that indeed for any e > 



and therefore 



lim lim F t ' e '^(\J[0,t]/t-j\ >e) = 

tiO tt+oo 



lim lim Ut e n = $i ■ 

40 tt+oo ' 



Now we estimate the entropy, arguing as in (3.23) 



, i — 1„/ , n i — 2(1 + 77) H(#- 
hm lim - H(za € „ a t ) < lim — — -— 

ej,0 tt+oo t ' " 40 TT + (p) + 7T~(p) 



-) + H(# H 



Now, recalling the definition of p in (3.17|)-(3.18), we have 



drk 



dn 



£—fa*-V*-t>i.)4 



so that 

H(tt" 



d(f> + d<t>~ d<t>+ (f)-(p) Pl 
log 4>~ (p) + / log p dir~ + log 



Pr 

Pl 



(Pr - Pi 



p 



TT (dp) . 
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As before, when e|0 we have — log 4> (p) + / log p drr 
e|0 



log a + log a = 0. Now, as 



P 



TT (dp) 



2 



(dp) 



Pr 



so that 



lmirlffr \<h" 

40 

and analogously 



H( 



■»^'l)-^^(l'-« fe -^ 



limH(r 
40 



log 



!h 



Therefore, by (3.25), 



limlim lim - H(vt e n I Pt) < Iolk 
T74.O ej.0 tf+oo t 



(Pr-Pl)(T l -T r ) 



(J3l-Pr)T l . 



(P R -h)\j\=I(j). 



Then again, there exists a map t h- y (e(t),r/(t)) vanishing as t f +oo such that vi := 
u t,e(t),r)(t) ~* S i and lim * t_1 H( ^' I /•**) - J ^)- D 

4. Scaling limit 
It is natural to consider the difference of temperatures as a variable and we introduce: 



1 



1 



H\r,T):=f[X, T+VT T ,. 



recall (3.7), and the corresponding Legendre transform: 



Then we have 



l(j,r,T) = sup{Aj-^(A, r,T)}. 

A 



e~ 2 l(ej, st, T) = sup{Aj - e~ l J : (eX, er, T)}. 

A 



The central result necessary to study the scaling limit is the following 



tU 



Proposition 4.1. Let K = (^) 2 - 
(1) J/r^O, then 



Xkt + kX 2 T 2 if At > 0, 



]im£- 2 F(e\,£T,T)=H(\,T,T):= I if AtG[-t 2 ,0], 

1 -(\ + t)kt + k(t + \) 2 T 2 if \t<-t 2 



(4.1) 



and 



lime- 2 X( £ j, £ r,T) = G(j,T,T) := { 



■J- K >r -r - . 2 

4kT 2 1X 3 T > „ KT 



f (j-«r) a 

4kT 2 x± J' ' 

if j'r G [0, kt 2 ] 
II if jrG[ _ Kr 2 )0] 



2T 2 



(4.2) 



£$£ if ir < -kt 2 . 



(2) Ifr = 0, lim £40 e" 2 ^(eA, er,T) = kA 2 T 2 and lim £40 £~ 2 ^J,£T,T) = ^. 
77ie above convergences are uniform in (X,t,T), resp. (j,T,T). 

The plot of Q was given in figure 1 and we give a plot of % in Figure 6. 
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Figure 6: Plot of H as a function of A for k = r = T 2 = 1 



Proof. We assume that r > 0, the case r < being completely analogous. We show first 
that VA > 0, 



lwi£~ 2 F(e\,£T,T) 

We set g(e) := F(eX,et,T) and use the notation: 

g(0+) = ]im g(t), g(0~ 



-Xkt + kX 2 T 2 . 



limq(t). 
tto 



By proposition (3.2) and the implicit function theorem J-(sX,£T,T) is a smooth function 



of s for e > 0. Thus, we write for e > small enough, 

g(e) = g(0 + ) + e</(0+) + y</'(0+) + 0(e 3 ) 

uniformly in (A, r, T). Since .F(-,0, T) is continuous and .^(O, 0,T) = by proposition |3. 2 
we get lim^o -F(eA, et, T) = 0, i.e. g(0 + ) = 



(4.3) 



We compute now the derivatives g'(0 + ) and </'(0 + ) using expression (3.9) and the relation 
(3.11). In order to simplify notations we introduce the distribution over M^_: 



ip(vi,v 2 ) = PlPrViv 2 exp 



(/3L-A)f-(/3 fl + A)f 

v\ v 2 2 2 



>1 



V 



(4.4) 



where rj is chosen such that the distribution is normalized, i.e ( |3.9| ) holds. We denote by 
^p L ,p R ,\(h) expectation of a Borel function h with respect to ip. Below ry is chosen so that 
(3.9) holds with the rescaled variables, namely r} = rf = J-(eX,£T,T). 



Let us start with the first derivative: 

</(0+) = lim — F(eX, et, T) = lim ( — (eA, er, T)X + ^-(eA, er, T)r 



(4.5) 



Taking first the derivative with respect to A, we get from (3.9), 

dF dF dF 
dr\ dX dX 
Next one computes explicitly, 



0. 



dF 
3A 



1. 



(eA, 77, /3l(et), /3k (et)) = -E / 3 L(£T)i/3fl(eT)i£A (r; 2 - Vl ). 



(4.6) 
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with P l (t) = {T + t)- 1 and /3 r (t) = {T - r)" 1 . Since E Mfi { 



0, we get 



dF 

lim ^r( eA ' r ?'^( £T )'^( er )) = 
£4.0 dX 



by dominated convergence. Taking next the derivative with respect to r, we obtain from 

(4.7) 



(3.9), 



dF dF dF dp L dF dfi R n 



dr\ dr dfii dr df3 R dr 



and 



Similarly, 



, dF 
hm— — 
£40 df3 L 



, 9F 
hm — — 
£^o <9/3r 



(eA,r?,^ L (er),/3 jR (er)) 



(eX,7],(3 L (£T),(3 R (eT)) 



T-^E Mi0 (vl) 



r--E A/3 ,oK) 



0. 



0. 



(4.8) 
(4.9) 



Using those two expressions in ( |4.7| ) , we get 

9 



lim T^-HeA, er T) = 0. 
£4.0 ar 



Combining this with (4.6) in (4.5), we finally obtain that </(0 + ) = 0. 



We next compute the second derivative with respect to e 



</'(0 + ) 



d 2 
lim -— tF{e\, £t, T) 
£|0 de 2 



The first two derivatives in (4.10) have been computed in [18J, the result is: 

'd 2 T, 



lim I a\9 

£40 V d\ l 



(e\,eT,T)\ z + 



d 2 T 



eX, et, T)Xt = 2(kT 2 X 2 + kXt). 



(4.11) 



dXdr 
We show now that the second derivative with respect to r vanishes when e goes to zero. In 



(3.11), we write (3l and /3r as function of r and take derivatives with respect to r. 

2 



d 2 F 
dr 2 



d 2 F ( dF 



dr/ 2 \d(3ft 



-P 



h' 



dF 

Wl 



01 



dF 

dr] 



, , dF «3 , dF S 



+ 






5?? / 



-i 



9 2 F 2 9 2 F 2 

PL/5R - 1^2 PL ~ W^2 PR 



m 



dp 



R 



dF 
dr] 



(4.12) 



By (4.8) and (4.9) the first derivatives of F with respect to (3l and /3r vanishes when e \. 0. 



We compute now the second derivatives with respect to (3l and (3 R and show that each also 
vanishes when e \, 0, for instance, 

d 2 F 



lim '„ 
£40 dpi 



(eX,r],p L (£T),f3 R (ET)) 



(-T 2 + |E y9l/ 3,o(«i) - \tE PAq {v\)\ = 0. 



All other derivatives may be dealt with in the same way and thus 

d 2 F , 



lim ., 
£40 dr 2 



(eX,er,T) = 0. 
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Combining this with (4.11) in (4.10), we finally get for the second derivative, 

g"{0 + ) = 2{ K T 2 \ 2 - K \T). 



Plugging this last expression in d4.3h, we finally obtain the result (4.1) for A > 0. 



Let us now consider the case A < — r < 0. By the Gallavotti-Cohen symmetry relation 



(3.10) and the definition of J-(X,t,T), we obtain that for A < — r, 

£~ 2 F(eX, er, T) = e' 2 F{e(j - A), er, T), 



and therefore by (4.1) valid for A > 0, we get: 



lime-^(eA,er,r) = (X- t)kt + k{t - X) 2 T Z if A < -r. 
£4,0 



Finally, if A G [— r, 0], then F(eX, er, T) = by point (2) of Proposition [3^ Formula ( |47T[ ) 
is now completely proved. The case r = is exactly analogous. 



We want to prove now (4.2). We assume again that r > and we first note that for any 
e > 0, 

if j € [0, kt] 

2T 2 

This follows from Lemma V 



e" 2 X(ej, er, T) = sup{Aj - e~ 2 .F(ej, er, T)} 

A 



=& if j € [-KT, 0] 



(4.13) 



When j > «;r, one notices that the maximizer X(e,j) of e 2 Aj — J-"(eA,er, T) is such that 
< X(e,j) < +oo and is the solution of the implicit equation in the unknown A: 



£J 



dX 



(eX,£T,T). 



(4.14) 



By the implicit function theorem, X(s,j) is therefore a smooth function of e > 0. By 
performing computations similar to the above ones, namely expanding (4.14) in e one can 
show that A(0 + , j) = ^~^2 which is the maximizer of the expression Xj — %{X, r, T), where 



T~L has been defined in (4.1). For each e > 0, e J-{eX,£t,T) is a convex function (as a 



function of A), thus the convergence to %(A, r, T) is uniform and 

.-2 



Therefore, 



lime- 2 T(eX(e,j),er,T)=n(X(0 + ,j),r,T). 
e4-0 



lime- 2 l(ej, er, T) = A(0+, j)j - H(A(0+, j),r, T) = G(j,r,T). 
e4-0 



(4.15) 



The case j < —kt is obtained from the case j > kt by using the Gallavotti-Cohen symmetry 
of Lemma 13.31 □ 

5. Appendix: The generator 

To avoid too heavy notations, we consider the simpler case of a particle moving in the 
interval [0, 1[ with a positive velocity. When the particle reaches 1, it is absorbed and re- 
emitted in with a random positive velocity distributed with a density ip. This dynamics 
has been introduced in section 2 of |18| and we follow the notations introduced there. The 



statement of proposition 2.1 is a trivial adaptation of proposition 5.1 below. We want 
to compute the infinitesimal generator, or more precisely the infinitesimal action of the 
dynamics on a smooth function (which is not necessarily in the domain of the generator). 
In other words, our aim is to prove the following result, where we denote the law of r, by 

"0((ir) and the law of v i = l/r« by ip{du). 
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Proposition 5.1. For all f,g : [0,1] x R + 

derivatives: 



bounded with bounded continuous first 



d_ 
dt 



I dq / dpg(q,p)P t f(q,p) 

JR+ 



t=0 



dp dqg(q,p)pf q (q,p)+ dppg(l,p) <p(du) (/(0,u) - /(!,£>)) 



Proof. The law of t\ + ■ ■ ■ + r n is denoted as usual by the n-fold convolution ip n * and we 
recall that S n = Sq + T\ + • ■ ■ + r n . Then we can write 



Ptf(qo,Po) = 

OO y 

= l(t<5 )/(*)+PO*.l^) + 1 (t>So) X^f^Sn-i^S™)/ 

n=l ^ 

= l(t<5 )/(90+Poi,Po) + 

oo „ „ 

+ l( t >5o)E/ r {n - l) {ds) ^{dr)f 

n=1 J[0,t-S Q ] J]t-So-s,+oo[ 



t — S n -\ 1^ 

7"n 7~n 



t — ^o — S 1 



r r 

t — 5o — s 1 

r r 



= l(t<5o)/(9o+Po*,Po) + l(t>5 ) / U{ds) I ip{dr)f 

J[0,t-S } J]t-S -s,+oc[ 

where we set ^*°(ds) = <5o(<i.s) and 

oo „5 oo „5 

n=l J a n=\ ■' a 

The renewal measure U{ds) gives the average number of collisions in the time interval ds. 
We define accordingly 



h 



(t) ■= / dpdqg(q,p)t {t<So ^ p)) P t f(q,p) 

J[0,l[xR+ 

rl-tp 

dp / dqg(q,p)f(q + tp,p), 

Jo 



h{t):= I dpdqg(q,p)t( t >s ( qiP)) P t f(q,p) 

'[0,l[xR+ 

1-9 



dpdqg{q,p)t {q y 1 _ tp) 

[0,l[xR + JO 



t-- 



+oo 



U{ds) / 1>{dv)f 



p l 



+oo rv p rl 

ip{dv) / U(ds) dp dqg(q,p)f 

(i JO JR+ JlA(l-tp+sp) 



V ' V 



Let us take the derivative in t 



W) = Jt h(t) 



dpp 



l-tp 



dqg(q,p) f g (q + tp,p) - g(l - tp,p) f(l,p) 
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■ x) /•« r r 1 l ( t-^— 1 - s \ s 

iP{dv) / U(cfe) / dp / dq-g(q,p)f q ? ,- 

(I JO Jr+ JlA(l-tp+sp) v \ V V 



-oo nv 

+ 1 ip(dv) U(ds) I dppla-tp +8 p<i)g(l-tp + sp,p)f[0, 

IQ JO -/IB ■ \ v 



+ 



Let us let t — > 0+: 



A(o) 



dpp 



dqg(q,p) fq(q,p) ~ s(l,p) /(l,p) 



and since U{ds) = 5o(ds) + l]o,+oo[( s ) U(ds) 



h 

Therefore 



(0)=/ dppg(l,p) [ iP(dv)f(0,v- 1 )= f dppg(l,p) [ <p(du) f(0,u). 

Jr + Jr + Jr + Jr + 



d 
dt 



JR+ 



dpg(q,p)P t f(q,p) 



t=o 



dp dqg(q,p)pf q (q,p)+ dppg(l,p) <p(du) (/(0,u) - /(l,p)) 

Jo Jr+ Jr+ 



□ 



The infinitesimal action of the dynamics on a smooth function / : [0, 1] x M + \— > M. can 
therefore be written 

Lf(q,p)=P^f+ P S 1 (dq) [ <p(du)(f(0,v)-f(l,p)). 

cq Jr+ 



5.1. The formal adjoint. Let us suppose now that ip(dv) = ifr(v)dv, so that the law of 

Vi = 1/n is 

tp(du) = ip(u) du = 5 — du. 



Then we can rewrite the result of Proposition 5.1 as follows 



d 
dt 



I dq I dpg(q,p)P t f(q,p) 

o Jr + 



t=o 



dp 

-!-- / dp 



dqg q (q,p)pf(q,p)+p(g(l,p)f(hp) ~ s(0,p)/(0,p)) 



^(p)/(0,p) / duug{l,u)- dppf(l,p)g(l,p) 

Jr+ Jr+ 



dp / dqg q (q,p)pf(q,p) 
Jo 



dpf(0,p) 



pg(0,p)-<p(p) I duug(l,u) 

+ 
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and obtain an expression for the formal adjoint of L 

L*g(q,p) = -p^--6 (dq) pg(0,p) - <p(p) / duug(l,u) . 
dq L Jr+ 

A solution of the Fokker-Planck equation associated with the process (qt,Pt) must then 
satisfy the boundary condition 



P9(0,p) = <p(p) duug(l,u), Vp>0. 

JWL+ 

We can check that g(q,p) = 1 [0,11(9) i p(p)/(P'P) is a probability density solving the equation 
L*g = 0. On the other hand, an invariant measure must satisfy L*g = 0. Since L*g is 
a sum of two mutually singular measures, they must both vanish. Then g(q,p) = g{p) is 
constant in q and 

Pg{p) = { f{p) I duug(u), 



i.e. 

I X 1 <P(P) 

9(P) = 7 • 

Z p 
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